Scaling theory for the free-energy barrier to homogeneous nucleation of a non-critical 
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Homogeneous nucleation of a new phase near an Ising-like critical point of another phase tran- 
sition is studied. A scaling analysis shows that the free energy barrier to nucleation contains a 
singular term with the same scaling as the order parameter associated with the critical point. The 
total magnetisation of the nucleus scales as the response function and so it diverges. Vapour-liquid 
critical points are in the Ising universality class and so our results imply that near such a critical 
point the number of molecules in a nucleus of a another phase, such as a crystalline phase, diverges 
as the isothermal compressibility. The case where symmetry prevents coupling between the nucleus 
and the order parameter is also considered. 



I. INTRODUCTION 



It is rare for a condensed matter system to possess only 
one phase transition, which raises the question of how 
neighbouring phase transitions affect each other. There 
are many possible effects, and here we consider one of 
them, namely how nucleation associated with one phase 
transition is affected by the critical point of a second 
transition. Consider a first-order phase transition, tran- 
sition a, which occurs at a temperature T^. Above Tq, the 
high-temperature phase, call it anr, is the equilibrium 
phase while below it the low-temperature phase, a^T: is 
the equilibrium phase. Fig. |l| is a schematic of the tem- 
perature axis of the phase diagram. If we cool the high- 
temperature phase below Tq, then it ceases to be the equi- 
librium phase. At some degree of supercooling a nucleus 
of the low temperature phase will appear and grow into a 
macroscopic aLT phase. The initial step, the formation 
of a nucleus of the aLT phase, is an activated process — 
a free energy barrier must be surmounted 0. The rate 
at which nuclei overcome a barrier of height AF* scales 
as exp(— AF*/fcr) and so this rate is very sensitive to 
the barrier height Here we calculate AF* for nucle- 
ation near a critical point of another phase transition, 
which we call transition /3. Transition /? is Ising-like. 
So, we have two phase transitions: a strongly first order 
phase transition a, whose dynamics we are studying, and 
a transition with a critical point, transition f3, which lies 
at a temperature < Ta and so occurs within the uht 
phase when it is not the true equilibrium phase. At an 
Ising critical point the thermodynamic functions contain 
singular power-law terms with exponents which depend 
only on dimensionality |3[^. Below, we show that AF* 
also contains a singular power-law term with exponents 
which depend only on dimensionality. The singular part 
of AF* has the same form as that in the order parameter, 
so for example, it scales with temperature as \T — Tfj\^ , 
where (3 is the usual critical exponent ||, Q . This singu- 



lar part means that the derivative of AF* with respect to 
temperature diverges as the critical point is approached: 
the barrier to nucleation is very sensitive to changes in 
temperature near the critical point. Which is true for 
any system with an Ising-like critical point. 

The scaling analysis here follows on from mean-field 
calculations performed by the author ||,||, sec also 
Ref. §. Which in turn followed on from the pioneering 
simulations of ten Wolde and Frenkel [|j who were the 
first to appreciate that critical fluctuations could reduce 
AF* near the critical point of a metastable fluid-fiuid 
transition. The results of the mean-field calculations are 
consistent with the scaling analysis here; if we take the 
equations of the present work and set the exponents to 
their (incorrect) mean-field values we reproduce the ex- 
ponents found in Refs. [^^. Finally, see the work of 
Dixit and Zukoski |^ for an alternative approach to nu- 
cleation near a critical point; one where dynamics are 
explicitly but approximately included. 

Perhaps the easiest systems in which to observe the 
effects found here, are mixtures of liquids. There nucle- 
ation of either the vapour phase or a crystalline phase 
may occur near the continuation of a line of critical solu- 
tion temperatures, i.e., near to a critical point of liquid- 
liquid demixing. A line of critical solution temperatures 
does not stop when the liquid ceases to be the equilibrium 
phase but continues into the region where it is metastable 
with respect to a transition to a vapour or crystalline 
phase. It is this continuation that affects nucleation. 

The study of nucleation near a critical point was 
sparked by interest in the nucleation of the crystalline 
phase of globular proteins |lO|. At least some globu- 
lar proteins have a fluid-fiuid critical point within the 
fiuid-crystal coexistence region The fact that 

the predictions made here are universal is very useful 
with regard to protein crystallisation as the interactions 
between protein molecules are generally rather poorly un- 
derstood, particularly before they have been crystallised 
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and their structure determined [Q. We know that our 
predictions will apply to nucleation near any fluid-fluid 
critical point found in a protein solution, without needing 
to know anything about the protein-protein interactions. 



II. THEORY FOR NUCLEI WHICH COUPLE TO 
THE ORDER PARAMETER 

We will describe the critical point of transition /? us- 
ing the standard language of Ising spin systems, al- 
though of course our results also apply to critical points 
in fluids. The order parameter is m, the temperature 
t = (T — Tp)/Ti3, and the field is h and couples to the 
order parameter as —hm. The critical point is then at 
t = 0, h = 0, see Fig. ^ For fluid systems the external 
field h is the chemical potential minus that at coexis- 
tence, and TO is a density difference. We set the lattice 
spacing to a. For definiteness we set m to be positive in 
the nucleus. The symmetry between +m and — m means 
we can do so without loss of generality. 

The dynamics of a first-order phase transition starts 
with nucleation ^ , where a microscopic nucleus of the 
new phase first appears. The nucleus at the top of the 
barrier, with the highest free energy AF*, determines the 
rate and so wc will consider only this nucleus. It is gener- 
ally called the critical nucleus. Note that conventionally 
the word critical is used to denote both a continuous tran- 
sition and the nucleus at the top of the barrier although 
there is no connection between the two uses of the word. 
Although this terminology is a little unfortunate we will 
use it here and so we will be studying a critical nucleus 
near a critical point. 

For nucleation at a strongly first-order phase transition 
(which transition a is assumed to be), the nucleus must 
contain a compact core of spins near their state in the 
bulk a^T phase, as only in this state do they have a lower 
free energy than in the anx phase. For all strongly first- 
order transitions, the free energy will be higher for spins 
intermediate between the two phases, here anx and olt- 
It is this core which will grow to form the a^T phase. See 
for example Refs. which discuss the classical theory 
for nucleation which is based on this observation. 

Now, this core of spins will perturb the surrounding 
spins out to a distance of the correlation length ^. Thus, 
the volume perturbed is of order ^'^ and so diverges as 
the critical point is approached, d is the dimensional- 
ity. The volume around the core where to is small we 
term the fringe. As in our earlier mean-field calculations 
1^,^ we split the nucleus into a core and a fringe. This 
is illustrated schematically in Fig. ||. As a is a strongly 
first-order transition, the order parameter of transition f3, 
TO, in the core will be of order unity, i.e., not small, and 
therefore far above the values of to reached by critical 
fiuctuations. The critical fluctuations will however affect 
the fringe where to is small. We note that the compact 



core of spins resembles a small colloidal particle in that 
it is a compact, roughly spherical object, which interacts 
with and thus perturbs its surroundings. The problem 
of a colloidal particle in a near critical system has been 
considered, see Refs. |l§-|l7|- 



The core of the nucleus will be much smaller than the 
correlation length ^ near the critical point as it will be 
only a few lattice sites across whereas ^ diverges 
Thus, near the critical point, the core is a point-like per- 
turbation which couples to the order parameter to. It 
will couple to TO unless this is prevented by a symmetry. 
Point-like perturbations near a critical point are generic 
— all point-like perturbations to the order parameter 
TO perturb large length-scale critical fluctuations in the 
same way, differing only in a scale factor. Near the criti- 
cal point where ^ ^ a, scale invariance implies that any 
perturbation of size much less than ^ must cause a or- 
der parameter perturbation of the same form, because all 
such perturbations are simply interrelated by a change of 
length-scale. The presence of a nucleus at the origin 
must therefore produce (at large distances) the generic 
perturbation to to. We define the thermal average excess 
of TO(r) due to the presence of a nucleus at the origin as 



< TO„(r) >=< m{r) >„ - < m{r) >, 



(1) 



where < m{r) >„ (< m{r) >) is the thermal average 
value of m(r) when there is (is not) a nucleus at the ori- 
gin. Compare this to the 2-point correlation function, 
G(r), for the order parameter, defined by 

G(r) =< TO(0)TO(r) > - < to(0) >< m{r) > . (2) 

Note that G(r) is the thermal average of m{r) if m(0) is 
fixed to 1 . As a point perturbation is generic for the fluc- 
tuations on length scales ^ a, the function < TO„(r) > 
can only differ from G(r) by a scale factor when r ^ a. 
Calling this scale factor G, 



< TO„(r) GG(r) r > a. 



(3) 



G will depend on the properties of the core. This result 
has been found before not for a nucleus but for a small 
colloidal particle immersed in a near-critical fluid [p^-|l7| . 
Defining the total excess of m due to the presence of the 
critical nucleus as to*, we have 



TO 



/ 



dr < TOn(r) > 



(4) 



Splitting the integration at fc which satisfies a ^ Tc ^ ^, 
and using Eq. for the integrand of the integral for 
r > Tc we have 



dr < m„(r) > +C 



drG(r). (5) 
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Now, the response function x is equal to the integral over 
G(r) il 



j drG(r) = X, 



(6) 



and it contains a leading order singular part xs com- 
ing from fluctuations on a length-scale of the correlation 
length ^, and so coming from those parts of the integral 
where r = 0{£). Thus the second integral of Eq. (||) for 
m* contains xs'- Tn* is singular at the critical point. Near 
the critical point xs dominates x ^-nd so the term Cxs 
in the second integral of Eq. (|^) will dominate m* . Thus, 
we have for the leading order singularity in m*, mj. 



m*s = Cxs^C\t\-'>X±ih/\tf'), 



(7) 



where X± are scaling functions, and 7, /? and S are criti- 
cal exponents. is the scaling function above the tran- 
sition, t > 0, while X_ is the scaling function below the 
transition, t < 0. The second equality is obtained by 
substituting the standard scaUng expression for xs § • 
Close to the critical point, the total order parameter ex- 
cess associated with the nucleus, m* , diverges with the 
precisely the same scaling behaviour as the order param- 
eter response function, x- 

Having determined the size of the nucleus, we wish to 
know its free energy. We start by splitting the free en- 
ergy barrier into a nonsingular part A_Fjv5 plus a leading 
order singular part AF* , 



AF* = AF, 



NS 



AF^ 



(8) 



The leading order singular part, AFg, of AF* , is ex- 
pected to have the usual scaling form, of a power of |i| 
times a function of h/\t\^^ , 



AF*s = \t\^Y±{h/\t 



(9) 



with X an as yet unknown exponent. To determine x we 
note that the derivative of the free energy of the system 
F with respect to h is, 



(10) 



where M is the total order parameter m of the system. 
Thus, as AF* is the free energy with a nucleus at the 
origin minus that without a nucleus at the origin, then 
the derivative of A_F* with respect to h is simply the neg- 
ative of the total order parameter with a nucleus minus 
that without a nucleus lllOl, 



dAF* 
dh 



(11) 



This is for the complete AF* . The leading-order singu- 
lar term in AF*, denoted by AFj, will yield minus the 



leading-order singular term in m* 
with respect to h 

dAFt 



dh 



upon differentiation 



(12) 



The derivative with respect to h of scaling functions like 
Y±{h/\t\f^^), brings an additional \t\-'^^ factor. Thus for 
Eqs. (Q) and (|l^) to yield the correct exponent —7 for m* 
we require that x — [35 = —7. Using the relation between 
exponents /3i5 = /? -I- 7 [^,||] we obtain x = {3. So, setting 
X = (3 in Eq. (Q) we obtain 



AF* = \tfY±{h/\tf'), 



(13) 



the singular part of the free energy scales with the same 
exponent, /?, as the bulk order parameter. This is for 
Ising-like systems but we expect that analogous singu- 
lar terms will be found near critical points of transitions 
with order parameters which are not scalars, for example 
in the Heisenberg model where the order parameter is a 
three-dimensional vector. 

As /3 > 0, AFg = at the critical point, and so the 
universal singular part is zero. The free energy barrier 
AF* at the critical point is equal to some non-universal 
value. However, as < 1, the derivative of AF* with 
respect to temperature diverges as It]''"^ as the critical 
point is approached at /i = 0. The signs of the diver- 
gences as the critical point is approached from above and 
below are those of Y±{Q) which are unknown. However, 
on physical grounds we might expect critical fluctuations 
to reduce the barrier to nucleation, and this is what is 
found in the mean-field theory calculations for T > Tp 
Ref. Q, i.e., within mean-field theory it has been found 
that Y+{Q) > 0. If 1±(0) are both positive, then at fixed 
h = 0, the free energy barrier will have a local minimum 
aX T = Tfj. As the derivative of AF| diverges, near the 
critical point the variation in AF* is dominated by that 
in AFg, and so if AFg has a minimum at the critical 
point so does AF* . 

Along the critical isotherm, 

AF^ = I]sgn(;i)|/i| 
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t = 0, 



(14) 



where E is an amplitude and the scaling with h is fixed by 
the requirement that it be the same as that of m. From 
Eq. (|l^) we see that as m*g > then the h derivative 
of AF^ must be negative. For this to be true we must 
have that S < 0. Note that m*g > because we set m 
in the core to be > 0, if the nucleating phase has instead 
m < 0, then in that case S > 0. S is a function not only 
of the properties of the critical system but of the nucle- 
ating phase. At the critical point the rate of decrease 
of AF* (= TO*) diverges, so near the critical point the 
nucleation barrier is a rapidly decreasing function of h. 
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A. Critical amplitude ratios 

In the h = Q hmit, the scahng functions X± of xs, 
Eq. (0), simphfy 



xs = r±\t\ 



0, 



(15) 



where T± — -'^±(0). T± are critical amplitudes, and 
are non-universal. However, the critical amplitude ratio 
r+/r_ is universal 0,|o|. The fact that not only the 
exponents but the critical amplitude ratios are univer- 
sal means that we can even make a universal prediction 
about the ratio of the size of the nucleus just above the 
critical point to that just below. If we take the ratio 
™s(+l*l)/™s(^l^l); we see, using Eq. (|l|), that 



^i+\t\ 



(16) 



which is universal. As near the critical point dom- 
inates m*, we have that near the critical point: 1) the 
excess magnetisation associated with the nucleus is the 
same in each of the two coexisting phases at some t < 0, 
and 2) that the ratio of m* at \t\ to the common value 
of the excess m in the coexisting phases at —\t\ is uni- 
versal. In three dimensions r+/r_ = 4.9 |§,|2y] so the 
nucleus is approximately five times bigger \t\ above the 
transition than —\t\ below. Via the nucleation theorem 
pl|-p^ it is possible to estimate the size of the critical 
nucleus. Thus measurements on nucleation rates can po- 
tentially be used to estimate not only critical exponents 
but critical amplitude ratios. 



III. DYNAMIC SCALING FOR THE 
EQUILIBRATION TIME NEAR Tp 

Now, the singular part of the barrier, AFg, is uni- 
versal because it originates in long wavelength fluctu- 
ations, fluctuations with wavelengths ^ a, where the 
microscopic details of the system are irrelevant. How- 
ever, these long wavelength fluctuations are slow, firstly 
as a direct consequence of them having a long wavelength 
(this applies if m is a conserved variable Q), and secondly 
because as the critical point is approached the free energy 
becomes a very flat function of the order parameter which 
reduces the driving force which relaxes fluctuations. Dy- 
namic scaling theory yields for the characteristic relax- 
ation time Tn {—2^ over the characteristic frequency) of 
a small wavevector mode, k << a^^, near a critical point 

i,|2i 



TR = k-^D±ikth/\tf') 



(17) 



the relaxation time of a mode of wavevector k at temper- 
ature t and field h. The exponent z is a dynamic scaling 
exponent. Dynamic scaling exponents depend not only 



on the symmetry of the order parameter and the dimen- 
sionality as the static exponents, /3, 7, etc. do, but also 
on whether the order parameter m is conserved or not. If 
m is the magnetisation then it can vary locally without 
transport of m and so it is non-conserved, but if m is a 
density difference then as the total number of molecules 
is conserved it cannot vary in a volume of space without 
transport in or out of that volume, so it is then a con- 
served variable. In either case z is positive. Assuming 
that fiuctuations with wavelengths larger than the cor- 
relation length can be ignored we have that at a given 
t and h, the longest relevant relaxation time is that of 
a mode of wavevector k = 2tt/^. Calling the relaxation 
time of this mode we have, from Eq. (|l7|). 



=r(27r)-^i?±(2^,0) 



(18) 



where we have set h — for simplicity, the generalisation 
to non-zero fields is easy. The second line used the scal- 
ing for the correlation length ^ — (,±\t\~'^ for h = 0, and 
the amplitudes t± = C^(27r)~^i:)±(27r, 0). 

Our expression for the free energy barrier to nucle- 
ation, Eq. (H) , includes a singular part AFg coming from 
long-wavelength fiuctuations. This singular part is an 
equilibrium quantity: it is the free energy change of equi- 
librium fluctuations due to the nucleus. Note that here 
by equilibrium we mean that the fluctuations in m etc. 
have relaxed to equilibrium, the system is not at true 
equilibrium as it is in the anr phase below Tq,. Thus, 
our expression for the free energy barrier, Eq. (|8|), is only 
meaningful once a time much larger than has elapsed 
at that temperature and so the fluctuations in m have 
equilibrated. 

However, diverges as the critical point is approached 
so the nearer we are to the critical point the longer we 
have to wait before Eq. (^ gives the correct barrier to 
nucleation. For m conserved, z = A — rj and so putting in 
the values for 77 and v we have that scales as |t|~^-^: it 
diverges rather rapidly as the critical point is approached. 
For m not conserved, z is close to 2 and so putting in the 
values for 77 and f we have that scales as |i|~'^ '^: it 
still diverges but more slowly. So, the closer the critical 
point is approached the longer is the waiting time before 
our expressions for the nucleation barrier, and size of the 
nucleus are valid. This, together with the effect of grav- 
ity on long wavelength fluctuations will limit how closely 
the critical point can be approached and the behaviour 
predicted by Eqs. (0) and (^ observed. 



IV. NUCLEI WHICH DO NOT COUPLE TO THE 
ORDER PARAMETER AND NON-ISING 
UNIVERSALITY CLASSES 

Section |l| considered the case of a nucleus near an 
Ising-like transition, and in which the nucleus couples to 
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the order parameter. In this section we start by consid- 
ering nuclei near an Ising-like critical point which do not 
couple to the order parameter, and then move on to con- 
sider critical points in other universality classes. There 
are results for colloidal particles in fluids near Ising-like 
critical points where the particles do not couple to the 
order parameter This is so when the two phases 

which form below Tp are related by symmetry and the 
nucleus does not discriminate between these two phases. 
Now, whether or not a nucleus (or colloidal particle) cou- 
ples to the order parameter, it will in general couple to 
the energy density u. The perturbation to the energy 
density around a particle/nucleus is given by the energy 
correlation function. A scaling argument analogous to 
that in the previous section can be constructed to yield 
the result that the coupling to the energy density yields 
a singular part in the free-energy barrier which scales as 
the energy density u, i.e., 

AF^ = \t\''^Z±{h/\tf'), (19) 

where a is the specific-heat exponent, and Z± are scaling 
functions. The right-hand side of Eq. (^9|) is the appro- 
priate scaling for the energy density. Now, 1 — a > /?, 
and so if a nucleus couples to both m and u then the 
singular term from the m coupling, Eq. (p^, will domi- 
nate and we can neglect the singularity of Eq. (|l9|): we 
did this in the previous section. In general, we always 
expect the nucleus to couple to the energy density but 
the resulting singularity, Eq. ( |l9|) will be dominated by 
that from the m coupling unless the latter is absent in 
which case Eq. ([l9| ) should give the dominant singularity 
in the free energy barrier to nucleation. Note that the 
temperature derivative of the nucleation barrier AF* /T. 
varies as the heat capacity, i.e., as |t|~" at /i = 0. 

So, for a nucleus which does not couple to the order 
parameter, we still have a singularity in the free-energy 
barrier, Eq. (p^, but it is weak. In three dimensions, 
a = 0.11 and so the temperature derivative of the free- 
energy barrier diverges only as at h = 0: this 
weak divergence may be difficult to observe. Also, at 
h = 0, ii the nucleus does not couple to m, then m* — 
at /i = by symmetry. Thus, Eq. (pT|), the field deriva- 
tive of dAF*/dh = at /i = 0. 

As mentioned in the introduction, in fluid systems like 
protein solutions and liquid mixtures, there is no symme- 
try available to prevent the coupling of the nucleus to m. 
So the results of the previous section apply there, and the 
coupling to the energy density is subdominant. However, 
let us consider systems in universality classes other than 
that of the Ising model. We speculate that our finding 
of singular terms in the free-energy barrier to nucleation 
near a critical point is not restricted to systems with a 
scalar order parameter. Even if the order parameter is, 
say a two-dimensional unit vector as it is in the univer- 
sality class of the XY- model [||,^, then the nucleus of a 



new phase should couple to the critical fluctuations. For 
example, the fluid to superfluid transition in liquid "^He 
is in the XY universality class. Now, the vector order pa- 
rameter is the phase of a wavefunction and a nucleus of 
another phase, e.g., the vapour phase, will not couple to 
the phase of a quantum mechanical wavefunction: AF* 
does not depend on the local phase of the wavefunction. 
However, there is no reason for the nucleus not to couple 
to energy fluctuations near the fluid to superfluid transi- 
tion. If it does so then AF* will contain a singular term 
of the form of Eq. (|l|). 

The fluid to superfluid transition crosses the liquid- 
vapour transition Q|, and so nucleation of the vapour 
phase can occur near the continuation of the superfluid- 
fluid transition into conditions of pressure and tempera- 
ture such that the vapour phase is the equilibrium phase. 
However, for the XY model in three dimensions, a is 
(very small but) negative which means the heat capac- 
ity, and thus the temperature derivative of A_F*, has a 
cusp only: it does not diverge. This may make proving 
or disproving our speculation difficult. 

V. CONCLUSION 

We have considered the effect of a critical point on the 
nucleation of a new phase of another phase transition. 
The order parameter associated with the critical point is 
a scalar, putting the transition in the universality class of 
the Ising model. The critical point results in power-law 
singularities in both the free energy barrier to nucleation, 
AF* , and the size of the nucleus, m*. The singular term 
in AF* has the same scaling as the order parameter, and 
that in m* has the scaling of the response function of 
the order parameter. Thus, for example, as the critical 
point is approached at zero field, the singular term in 
AF* varies as \t\^ , where /3 is the usual critical exponent 
and t is the temperature minus that at the critical point. 
As /3 > 0, the singular term in the free energy barrier 
tends to zero as the critical point is approached. How- 
ever, derivatives of AF*, both with respect to the tem- 
perature and to the external field, diverge as the critical 
point is approached: the free-energy barrier to nucleation 
varies rapidly near the critical point. These predictions 
are universal, they apply to nucleation in any system in 
which nucleation of a new phase occurs near a critical 
point in the universality class of the Ising model. 

All the above holds given only that the core of the nu- 
cleus, where the nucleus resembles the new bulk phase 
that is nucleating, couples to the order parameter. It 
will do so unless prevented by a symmetry. This symme- 
try is clearly absent in fluid systems, in which the order 
parameter is a density and/or composition difference. In 
the previous section we speculated that there are also 
singular terms in AF* near critical points in other uni- 
versality classes. In other universality classes it is clear 
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that symmetry can prevent coupling to the order param- 
eter, in which case the singularity in AF* will be weaker. 

The relaxation time diverges as the critical point is ap- 
proached, see section [II. Thus, in experiment there is a 
limit to how close the critical point can be approached be- 
fore the relaxation time becomes prohibitively long. This 
may limit the minimum value of |t| that is achievable or 
it may be limited by some other factor. For example, 
the protein lysozyme has a fluid-fluid transition with a 
critical point within the fluid-crystal coexistence region 
and thus is a candidate for studying nucleation of 
a non-critical phase, here a crystalline phase, near an 
Ising-like critical point. Far from the critical point the 
characteristic relaxation time can be estimated as the 
time a lysozyme molecules takes to diffuse its own diam- 
eter. The single particle diffusion constant of lysozyme 
D = O(10~^°m^s) and its diameter is approximately 
4nm This gives a relaxation time of order 10~^s. Us- 
ing this as the order of magnitude of the amplitudes r±, 
we have that at h — 0, the relaxation time is of order 
10~^|t|~^'^s at a temperature t. The exponent of t is zi', 
with z taking its value for a conserved order parameter. 
Now, assuming careful temperature control the critical 
temperature can be approached to within, say, 0.1°C. As 
the effective interactions between the molecules in solu- 
tion vary strongly with temperature 0.1°C should corre- 
spond to \t\ = 0(10-2 - 10-3), not \t\ = 4x 10-^ as 
it would for temperature independent interactions. At 
\t\ = 10-3, we have a relaxation time — 0(ls): still 
rather short. Thus, in experiments on globular protein 
molecules like lysozyme, limitations on how well the tem- 
perature can be controlled are the limiting factor; the re- 
laxation time is always well within the range accessible 
in experiment. 

It is a pleasure to acknowledge discussions with R. 
Evans, A. Parry and P. Upton. This work was supported 
by EPSRC (GR/N36981). 
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FIG. 1. Schematic phase behaviour along the temperature axis, at /i = 0. The shaded (hatched) region denotes the 
temperature range over which the high (low) temperature phase of transition a is the equilibrium phase. The temperatures of 
the two transitions are labelled by and T^. 




FIG. 2. Schematic of a nucleus of the ordered phase of transition a near transition /3. The core of the ordered phase of 
transition a is solid black, and the perturbation this causes iu the surroundings is the shaded circle of radius the correlation 
length ^. The spherical boundary with radius Tc is denoted by the dashed circle. 
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